Choi ef al. Breast Cancer Research 201 3, 1 5:R78 
http://breast-cancer-research.eom/content/15/5/R78 



Breast Cancer 

RESEARCH 



RESEARCH ARTICLE 



Open Access 



Metabolic interaction between cancer cells and 
stromal cells according to breast cancer 
molecular subtype 



Junjeong Choi 1 , Do Hee Kim 2 , Woo Hee Jung 2 and Ja Seung Koo 



2* 



Abstract 

Introduction: The aim of this study was to investigate the differential expression of markers related to metabolic, 
mitochondrial and autophagy status in different molecular subtypes of breast cancer. 

Methods: Using tissue microarray sections generated from 740 cases of breast cancer, we performed 
immunohistochemical staining for Glut-1, CAIX, MCT4, ATP synthase, glutaminase, BNIP3, Beclin-1, LC3A, LC3B and 
p62. Based on the immunohistochemical expression of estrogen receptor (ER), progesterone (PR), HER2, and Ki-67 
labeling index, the cases were classified into luminal A, luminal B, HER2 and triple-negative breast cancer (TNBC). 
We further classified metabolic phenotypes of tumors according to glycolytic status by assessing Glut-1 and CAIX 
expression as follows: Warburg type: tumor (glycolysis type), stroma (nonglycolysis type); reverse Warburg type: 
tumor (nonglycolysis type), stroma (glycolysis type); mixed type: tumor (glycolysis type), stroma (glycolysis type); 
and null type: tumor (nonglycolysis type), stroma (nonglycolysis type). 

Results: Expression of Glut-1, MCT4 and LC3A was highest in TNBC and lowest in luminal A type (P < 0.001). 
Tumors were classified into 298 Warburg type (40.3%), 54 reverse Warburg type (7.3%), 62 mixed type (8.4%) and 
326 null type (44.0%). The mixed type had a higher histologic grade, ER negativity, PR negativity and Ki-67 index, 
whereas the null type showed lower histologic grade, ER positivity, PR positivity and Ki-67 index (P < 0.001). TNBC 
constituted the major portion of Warburg and mixed types, and luminal A consisted mainly of reverse Warburg and 
null types (P < 0.001). 

Conclusion: Breast cancer is heterogeneous in its metabolic status, and therefore it can be classified into various 
metabolic phenotypes. Specifically, the Warburg and mixed types had strong associations with TNBC, whereas 
reverse the Warburg and null types had associations with the luminal type, suggesting a correlation between 
metabolic phenotype and the biology of breast cancer. 



Introduction 

The metabolism of malignant tumors is generally explained 
by the Warburg effect theory, which describes the meta- 
bolic shift from mitochondrial oxidative phosphorylation 
(OXPHOS) to glycolysis in tumors [1]. Breast cancer is 
known to be heterogeneous, and the interaction be- 
tween tumor cells and adjacent stroma is expected to 
have significant roles in tumor growth and progression. 
This kind of complex interaction may also exist in the 
metabolic processes of the tumor. Previous studies 



suggest a unique metabolic interaction between tumor 
cells and the stroma of breast cancer, known as the reverse 
Warburg effect theory [2-5]. According to this theory, 
reactive oxygen species (ROS), such as nitric oxide 
(NO), generated by tumor cells bring oxidative stress 
to the stromal cells, leading to mitochondrial dysfunc- 
tion, autophagy (mitophagy) and increased aerobic gly- 
colysis through hypoxia-inducible factor la (HIF-la) 
and nuclear factor kB (NF-kB). Lactate generated by 
stromal cell glycolysis enters tumor cells and promotes 



* Correspondence: kjsl976@yuhs.ac 

department of Pathology, Yonsei University College of Medicine, 50 Yonsei-ro, 

Seodaemun-gu, Seoul 120-752, South Korea 

Full list of author information is available at the end of the article 



© 201 3 Choi et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BlOlVlGCl C^ntrBl Commons Attribution License (http//creativecomn ons.org/lii :nsi 'by/2.0), 1 hich perm ts unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Choi et al. Breast Cancer Research 201 3, 1 5:R78 
http://breast-cancer-research.eom/content/15/5/R78 



Page 2 of 20 



tumor cell growth and survival through efficient gener- 
ation of ATP by OXPHOS in the mitochondria. Cancer- 
associated fibroblasts (CAFs), stromal cells with loss of 
caveolin-1 expression, have been implicated in this inter- 
action in breast cancer because loss of caveolin-1 results 
from increased proteolysis by autophagy [3,5-7]. In addition 
to the reverse Warburg effect theory, it has been reported 
that certain types of tumor cells generate ATP through 
glycolysis as well as OXPHOS, which suggests various 
features of tumor metabolism [8,9] . Proteins involved in 
metabolism, mitochondrial function and autophagy may 
be differentially expressed in both tumor and stromal cells 
according to the aforementioned theories. These differences 
are summarized in Table 1. 

Because breast cancer is heterogeneous with respect to 
clinical, histopathological and molecular features, several 
subclassifications have been investigated to stratify tumors 
with similar characteristics. Gene expression profiles have 
enabled molecular classification of tumors into luminal A, 
luminal B, HER2, normal breast-like and basal-like types 
[10-12]. Studies have also revealed differences in histo- 
logical and clinical manifestation between different molecu- 
lar subtypes, such as therapeutic response and prognosis. 
We hypothesized that the metabolic interaction between 
tumor cells and stroma may differ according to the mo- 
lecular subtypes of breast cancer. Because there are lim- 
ited studies regarding this question, the aim of our present 
study was designed to investigate the differential expression 
of markers for metabolic, mitochondrial and autophagy 
status in different molecular subtypes of breast cancer. 

Methods 

Patient selection 

Patients diagnosed with invasive breast cancer treated by 
surgical resection during the period from January 2002 
to December 2006 were included in this study. Patients 
who received preoperative neoadjuvant chemotherapy or 
hormonal treatment were excluded. This study was ap- 
proved by the Institutional Review Board (IRB) of Yonsei 
University Severance Hospital. The IRB exempted the 
informed consent from patients. A breast pathologist 
(JSK) retrospectively reviewed the histology of all cases 
using hematoxylin and eosin (H&E)-stained slides. The 
histological grade was assessed using the Nottingham 



grading system [13]. Clinicopathologic parameters eval- 
uated in each case included patient age at initial diagnosis, 
lymph node metastasis, tumor recurrence, distant metasta- 
sis and patient survival. 

Tissue microarray 

On H&E-stained slides of tumors, a representative area was 
selected and the corresponding spot was marked on the 
surface of the paraffin block. Using a biopsy needle, the se- 
lected area was punched out and a 3-mm tissue core was 
placed into a 6 x 5 recipient block. Tissue from the invasive 
tumor was then extracted. More than two tissue cores were 
extracted to minimize extraction bias. Each tissue core 
was assigned a unique tissue microarray (TMA) location 
number that was linked to a database containing other 
clinicopathologic data. 

Immunohistochemistry 

The antibodies used for immunohistochemistry (IHC) in 
this study are shown in Table 2. Formalin-fixed, paraffin- 
embedded (FFPE) tissue sections from the TMA were 
prepared for IHC. Briefly, 5-um-thick sections were 
obtained using a microtome, transferred into adhesive 
slides and dried at 62°C for 30 min. After incubation 
with primary antibodies, immunodetection was performed 
with biotinylated anti-mouse immunoglobulin, followed by 
peroxidase-labeled streptavidin using a labeled streptavidin 
biotin kit with 3,3 -diaminobenzidine chromogen as 
the substrate. The primary antibody incubation step was 
omitted in the negative control. A positive control was in- 
cluded for each experiment: glucose transporter 1 (Glut-1): 
esophageal carcinoma; carbonic anhydrase IX (CAIX): renal 
carcinoma; monocarboxylate transporter 4 (MCT4): BCL2/ 
adenovirus E1B 19-kDa interacting protein 3 (BNIP3), 
kidney tissue; Beclin-1: breast tissue; microtubule-associated 
protein 1 light chain 3a (LC3A): brain tissue; microtubule- 
associated protein 1 light chain 3(3 (LC3B): brain tissue, 
p62: spleen tissue; ATP synthase: heart tissue; and glu- 
taminase: liver tissue. Slides were counterstained with 
Harris hematoxylin. 

Interpretation of immunohistochemical staining 

All immunohistochemical markers were assessed by light 
microscopy. Pathologic parameters such as ER, PR and 



Table 1 Comparison of metabolism type, mitochondrial status and autophagy status between the Warburg effect 
theory and reverse Warburg effect theory 3 





Warburg effect theory 
Cancer cell 


Reverse Warburg effect theory 
Cancer cell Stromal cell 


Metabolism 


Glycolysis 


OXPHOS 


Glycolysis 


Mitochondrial status 


Dysfunctional 


Functional 


Dysfunctional 


Autophagy status 


Not included 


Not activated 


Activated 



'OXPHOS, oxidative phosphorylation. 
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Table 2 Source, clone and dilution of antibodies used in this study 3 



Antibody 



Clone 



Dilution 



Company 



Molecular subtype-related 
ER 
PR 

HER2 
Ki-67 
Glycolysis-related 
Glut-1 
CAIX 
MCT4 

Mitochondrial status-related 
BNIP3 

Mitochondrial metabolism-related 

ATP synthase 

Glutaminase 
Autophagy-related 

Beclin-1 

LC3A 

LC3B 

p62 



SP1 
PgR 

Polyclonal 
MIB-1 

SPM498 

Polyclonal 

Polyclonal 

Ana40 

15H4C4 
Polyclona 

Polyclona 
EP1528Y 
Polyclona 
SQSTM1 



1:100 
1:50 
1:1,500 
1:150 

1:200 
1:100 
1:100 

1:100 

1:100 
1:100 

1:100 
1:100 
1:100 
1:100 



Thermo Scientific, Waltham, MA, USA 
Dako Denmark AS, Glostrup, Denmark 
Dako Denmark AS, Glostrup, Denmark 
Dako Denmark AS, Glostrup, Denmark 

Abeam, Cambridge, UK 
Abeam, Cambridge, UK 

Santa Cruz Biotechnology, Santa Cruz, CA, USA 

Abeam, Cambridge, UK 

Abeam, Cambridge, UK 
Abeam, Cambridge, UK 

Abeam, Cambridge, UK 
Abeam, Cambridge, UK 
Abeam, Cambridge, UK 
Abeam, Cambridge, UK 



a BNIP3, BCL2/adenovirus E1B 19-kDa interacting protein 3; CAIX, carbonic anhydrase IX; ER, estrogen receptor; Glut-1, 
protein 1 light chain 3a; LC3B, microtubule-associated protein 1 light chain 3(5; MCT4, monocarboxylate transporter 4, 



glucose transporter 1; LC3A, microtubule-associated 
PR, progesterone receptor. 



HER2 expression were obtained from each patient's 
pathologic report. A cutoff value of 1% or more positively 
stained nuclei was used to define ER and PR positivity 
[14]. HER2 staining was analyzed according to the American 
Society of Clinical Oncology (ASCO)/College of American 
Pathologists (CAP) guidelines using the following cat- 
egories: 0 = no immunostaining; 1+ = weak, incom- 
plete membranous staining, less than 10% of tumor 
cells; 2+ = complete membranous staining, either uniform 
or weak in at least 10% of tumor cells; and 3+ = uniform 
intense membranous staining in at least 30% of tumor cells 
[15]. HER2 immunostaining was considered positive when 
strong (3+) membranous staining was observed, whereas 
cases with 0 to 1+ were regarded as negative. Cases show- 
ing 2+ HER2 expression were evaluated for HER2 amplifi- 
cation by fluorescence in situ hybridization (FISH). 

Glut-1, CAIX, BNIP3, MCT4, Beclin-1, LC3A, LC3B 
and p62 immunohistochemical staining was evaluated 
on the basis of the proportion of stained cells and im- 
munostaining intensity. The proportion of stained cells 
was graded 0 (negative), 1 (less than 30% positive) or 2 
(more than 30% positive). Immunostaining intensity was 
graded as 0 (negative), 1 (weak), 2 (moderate) or 3 (strong). 
The scores for the proportion of stained cells and staining 
intensity were multiplied to provide a total score: negative 
(0 or 1) or positive (2 through 6). Ki-67 labeling indices 
(Lis) were scored by counting the number of positively 



stained nuclei and expressed as a percentage of total 
tumor cells. 

Fluorescence in situ hybridization analysis 

Before FISH analysis, invasive tumors were examined on 
H&E-stained slides. FISH was subsequently performed 
on the confirmed tumor. FISH was performed using the 
PathVysion HER-2 DNA Probe Kit (Abbott Molecular, 
Abbott Park, IL, USA) according to the manufacturer's 
instructions. HER2 gene copy number on the slides was 
evaluated using an epifluorescence microscope (Olympus, 
Tokyo, Japan). At least 60 tumor cell nuclei in three 
separate regions were investigated for HER2 and chromo- 
some 17 signals. HER2 gene amplification was determined 
according to the ASCO/CAP guidelines [15]. An absolute 
HER2 gene copy number lower than 4 or a HER2 gene/ 
chromosome 17 (chrl7) copy number ratio (HER2/chrl7 
ratio) less than 1.8 was considered HER2-negative. An 
absolute HER2 copy number between 4 and 6 or a HER2/ 
chrl7 ratio between 1.8 and 2.2 was considered HER2- 
equivocal. An absolute HER2 copy number greater than 
6 or a HER2/chrl7 ratio higher than 2.2 was considered 
//£7?2-positive. 

Tumor phenotype classification 

In this study, we classified breast cancer phenotypes 
according to the IHC results for ER, PR, HER2 and Ki-67 
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LI. FISH results for HER2 were as follows [16]: luminal 
A type: ER- and/or PR-positive, HER2-negative Ki-67 LI 
less than 14%; luminal B type: (HER2-negative) ER- and/or 
PR-positive, HER2-negative and Ki-67 LI greater than or 
equal to 14% and (HER2-positive) ER- and/or PR-positive 
and HER2 overexpressed and/or amplified; HER2 type: ER- 
and PR-negative and HER2 overexpressed and/or amplified; 
TNBC type: ER-, PR- and HER2-negative. 

Classification of tumor metabolic subtypes 

We also classified cases based on the results of immuno- 
histochemical staining for metabolism-related proteins 
as follows: glycolysis type: Glut-1- and/or CAIX-positive; 
nonglycolysis type: Glut-1- and CAIX-negative; dysfunc- 
tional mitochondrial type: BNIP3-positive [17,18]; func- 
tional mitochondrial type: BNIP3-negative; activated 
autophagy type: positive for two or more markers from 
among Beclin-1, LC3A, LC3B and p62; and nonactivated 
autophagy type: positive for less than two markers 
from among Beclin-1, LC3A, LC3B and p62. We further 
classified the metabolic phenotypes of breast cancer as 
follows: Warburg type: tumor (glycolysis type), stroma 
(nonglycolysis type); reverse Warburg type: tumor 
(nonglycolysis type), stroma (glycolysis type); mixed type: 
tumor (glycolysis type), stroma (glycolysis); and null type: 
tumor (nonglycolysis type), stroma (non-glycolysis type). 

Laser microdissection and protein extraction from 
formalin-fixed, paraffin-embedded tissues 

To acquire tumors and tumor stroma, laser microdissec- 
tion was performed with hematoxylin-stained, uncovered 
slides generated with FFPE blocks (LMD 6500; Leica, 
Wetzlar, Germany). Five cases per molecular subtype of 
breast cancer were microdissected. Protein extraction 
from microdissected FFPE tissues was performed using 
the Qproteome FFPE Tissue Kit (QIAGEN, Hilden, 
Germany). Briefly, microdissected FFPE tissues were 
deparaffinized in xylene and rehydrated in a graded 
series of alcohol. Afterward, the samples were mixed 
with FFPE extraction buffer EXB Plus (100 ul per sam- 
ple; QIAGEN), incubated at 100°C for 20 min, at 80°C 
for 2 h and then centrifuged for 15 min at 14,000 x g at 
4°C. The protein concentrations in the supernatant were 
determined using the Bradford assay (Bio-Rad Labora- 
tories, Hercules, CA, USA). 

Western blot analysis 

Total protein (20 ug) from each sample was mixed with 
Laemmli sample buffer and heated at 100°C for 5 min. 
It was then loaded into individual wells, resolved by 
8% SDS-PAGE and electroblotted onto nitrocellulose 
membranes (GE Healthcare Life Sciences, Pittsburgh, 
PA, USA). Membranes were blocked in 5% nonfat dry 
milk in Tris-buffered saline with Tween 20 (TBS-T), 



then incubated with antibodies to Glut-1, CAIX, ATP 
synthase, glutaminase, MCT-4, LC3A and p62 overnight 
at 4°C. The membranes were washed with TBS-T and 
then probed with peroxidase-conjugated goat anti-rabbit 
immunoglobulin G (1:2,000; Santa Cruz Biotechnology, 
Santa Cruz, CA, USA) for 1 h at room temperature. 
Washing was repeated and the membranes were developed 
with an enhanced chemiluminescence agent (Amersham/ 
GE Healthcare Life Sciences, Little Chalfont, UK). Band 
densities were measured using TINA image software 
(raytest, Straubenhardt, Germany). 

Statistical analyses 

Data were processed using SPSS for Windows version 
12.0 software (SPSS Inc, Chicago, IL, USA). Student's 
t-test and Fisher's exact test were used to examine any 
differences in continuous and categorical variables, respect- 
ively. Significance was assumed when P < 0.05. Kaplan- 
Meier survival curves and logrank statistics were employed 
to evaluate time to tumor recurrence and time to survival. 
Multivariate regression analysis was performed using the 
Cox proportional hazards model. 

Results 

Patients' characteristics according to tumor phenotype 

The clinicopathologic characteristics of 740 patients, 
comprising 298 (40.3%) cases of luminal A, 166 (22.4%) 
cases of luminal B, 69 (9.3%) cases of HER2 type and 
207 (28%) cases of TNBC type, are summarized in Table 3. 
TNBC type had the highest histologic grade, tumor 
stage and Ki-67 LI (P < 0.001, P = 0.002 and P < 0.001, 
respectively). In addition, HER2 and TNBC types had 
higher incidences of tumor recurrence and patient death 
than other types (P < 0.001). 

Expression of metabolism-related proteins according to 
tumor phenotype 

The differential expression of metabolism-related proteins 
according to breast cancer phenotype is summarized in 
Table 4. Tumor expression of Glut-1, MCT4 and LC3A 
was highest in TNBC and lowest in the luminal A type 
(P < 0.001). Stromal expression of CAIX and MCT4 
and tumor expression of cytoplasmic p62 was highest 
in HER2 type and lowest in luminal A type (P = 0.032, 
P < 0.001 and P < 0.001, respectively). Tumor expression 
of CAIX and LC3B was highest in TNBC and lowest in 
luminal B type (P = 0.008 and P = 0.013, respectively). 
HER2 type showed the highest tumor and stromal ATP 
synthase expression (P = 0.027 and P < 0.001, respectively) 
and stromal glutaminase expression (P = 0.001), whereas 
luminal A type showed the lowest expression of those 
markers. Expression of stromal LC3A and tumor expres- 
sion of nuclear p62 were highest in luminal A and lowest in 
TNBC {P < 0.001). 
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Table 3 Clinicopathologic characteristics of patients according to breast cancer phenotype 3 



Parameters 


Total 
(W = 740) (%) 


Luminal A 
(n = 298) (%) 


Luminal B 
(n = 166) (%) 


HER2 

(n = 69) (%) 


TNBC 

(n = 207) (%) 


P-value 


Age (years, mean ± SD) 


49.7 ±11.0 


50.6 ± 10.5 


48.5 ±10.1 


52.8 ± 9.8 


48.4 ± 1 2.4 


0.007 


Histologic grade 












<0.001 


I 


118 (15.9) 


90 (30.2) 


18 (10.8) 


1 (14) 


9 (4.3) 




II 


373 (50.4) 


1 80 (60.4) 


90 (54.2) 


35 (50.7) 


68 (32.9) 




III 


249 (33.6) 


28 (9.4) 


58 (34.9) 


33 (47.8) 


130 (62.8) 




Tumor stage 












0.002 


T1 


358 (48.4) 


166 (55.7) 


86 (51.8) 


31 (44.9) 


75 (36.2) 




T2 


367 (49.6) 


125 (41.9) 


78 (47.0) 


37 (53.6) 


127 (61.4) 




T3 


1 5 (2.0) 


7 (2.3) 


2 (1.2) 


1 (14) 


5 (2.4) 




Nodal stage 












0.041 


NO 


436 (58.9) 


1 68 (56.4) 


91 (54.8) 


42 (60.9) 


135 (65.2) 




Nl 


200 (27.0) 


90 (30.2) 


43 (25.9) 


13 (18.8) 


54 (26.1) 




N2 


66 (8.9) 


27 (9.1) 


17 (18.5) 


10 (14.5) 


12 (5.8) 




N3 


38 (5.1) 


1 3 (4.4) 


1 5 (9.0) 


4 (5.8) 


6 (2.9) 




Estrogen receptor status 












<0.001 


Negative 


286 (38.6) 


5 (1.7) 


5 (3.0) 


69 (100.0) 


207 (100.0) 




Positive 


454 (61.4) 


293 (98.3) 


161 (97.0) 


0 (0.0) 


0 (0.0) 




Progesterone receptor status 












<0.001 


Negative 


371 (50.1) 


50 (16.8) 


46 (27.7) 


69 (1 00.0) 


207 (100.0) 




Positive 


369 (49.9) 


248 (83.2) 


120 (72.3) 


0 (0.0) 


0 (0.0) 




HER2 status 












<0.001 


0 


290 (39.2) 


1 08 (36.2) 


23 (13.9) 


0 (0.0) 


159 (76.8) 




1 + 


186 (25.1) 


118 (39.6) 


33 (20.0) 


0 (0.0) 


35 (16.9) 




2+ 


142 (19.2) 


72 (24.2) 


41 (24.7) 


16 (23.2) 


13 (6.3) 




3+ 


122 (16.5) 


0 (0.0) 


69 (41 .6) 


53 (76.8) 


0 (0.0) 




Ki-67 LI (%, mean ± SD) 


18.1 ± 19.2 


4.7 ± 3.7 


19.7 ± 12.7 


19.5 ± 12.5 


35.6 ± 23.7 


<0.001 


Tumor recurrence 


69 (9.3) 


15 (5.0) 


12 (7.2) 


11 (15.9) 


31 (15.0) 


<0.001 


Patient's death 


67 (9.1) 


14 (4.7) 


11 (6.6) 


12 (17.4) 


30 (14.5) 


<0.001 


Duration of clinical follow-up (months, mean ± SD) 


70.2 ±31.7 


72.7 ± 30.0 


70.3 ± 30.3 


67.1 ± 35.8 


67.8 ± 33.8 


0.291 



a LI, labeling index; TNBC triple-negative breast cancer. 



Correlation between metabolism-related proteins and 
clinicopathologic factors 

The correlation between expression of metabolism-related 
proteins and clinicopathologic parameters is summarized 
in Table 5. Tumor expression of Glut-1 was associated 
with higher histologic grade (P < 0.001), ER negativity 
(P < 0.001), higher T stage (P < 0.001) and higher Ki- 
67 LI (P < 0.001), whereas CAIX was associated with 
higher Ki-67 LI {P < 0.001). Stromal ATP synthase expres- 
sion was associated with HER2 positivity (P < 0.001), 
and stromal glutaminase expression was associated with 
higher KI-67 LI (P = 0.021). Tumor expression of MCT4 
was associated with higher histologic grade {P < 0.001), 
ER negativity (P < 0.001), PR negativity (P < 0.001), higher 
T stage (P < 0.001) and higher Ki-67 LI (P < 0.001). 



Stromal expression of MCT4 was associated with higher 
histologic grade (P < 0.001), ER negativity (P < 0.001), 
PR negativity {P < 0.001), HER2 positivity (P < 0.001) 
and higher Ki-67 LI (P < 0.001). Tumor expression 
of LC3A was associated with higher histologic grade 
(P < 0.001), ER negativity (P < 0.001), PR negativity 
(P < 0.001), HER2 negativity (P < 0.001) and higher Ki-67 
LI (P < 0.001). In contrast, stromal expression of LC3A 
was associated with ER positivity (P < 0.001), PR posi- 
tivity (P < 0.001) and lower Ki-67 LI (P = 0.032). Tumor 
expression of cytoplasmic p62 was associated with 
HER2 positivity {P < 0.001), whereas nuclear p62 was 
associated with lower histologic grade (P < 0.001), ER 
positivity (P < 0.001), PR positivity (P < 0.001) and lower 
Ki-67 LI (P < 0.001). 
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Table 4 Expression of metabolism-related proteins according to breast cancer phenotype 3 



Parameters 



Total Luminal A Luminal B HER2 TNBC P-value 

(N = 740 ) (%) (n = 298)(%) (n=166)(%) (n = 69) (%) (n = 207) (%) 



Glut-1 in tumor 

Negative 504(68.1) 

Positive 236(31.9) 

Glut-1 in stroma 

Negative 724 (97.8) 

Positive 1 6 (2.2) 

CAIX in tumor 

Negative 520 (70.3) 

Positive 220 (29.7) 

CAIX in stroma 

Negative 627 (84.7) 

Positive 113(15.3) 

ATP synthase in tumor 

Negative 30(4.1) 

Positive 710(95.9) 

ATP synthase in stroma 

Negative 570 (77.0) 

Positive 1 70 (23.0) 

Glutaminase in tumor 

Negative 219(29.6) 

Positive 521 (70.4) 

Glutaminase in stroma 

Negative 495 (66.9) 

Positive 245 (33.1) 

BNIP3 in tumor 

Negative 504(68.1) 

Positive 236(31.9) 

BNIP3 in stroma 

Negative 700 (94.6) 

Positive 40 (5.4) 

MCT4 in tumor 

Negative 540 (73.0) 

Positive 200 (27.0) 

MCT4 in stroma 

Negative 418(56.5) 

Positive 322 (43.5) 

Cytoplasmic Beclin-1 

Negative 406 (54.9) 

Positive 334(45.1) 

Nuclear Beclin-1 

Negative 666 (90.0) 

Positive 74(10.0) 



260 (87.2) 
38 (12.8) 

296 (99.3) 
2 (0.7) 

217 (72.8) 
81 (27.2) 

264 (88.6) 
34 (1 1 .4) 

20 (6.7) 
278 (93.3) 

256 (85.9) 
42 (14.1) 

85 (28.5) 
213 (71.5) 

223 (74.8) 

75 (25.2) 

206 (69.1) 
92 (30.9) 

281 (94.3) 
17 (5.7) 

253 (84.9) 
45 (15.1) 

222 (74.5) 

76 (25.5) 

169 (56.7) 
129 (43.3) 

262 (87.9) 
36 (12.1) 



1 24 (74.7) 
42 (25.3) 

1 62 (97.6) 
4 (2.4) 

127 (76.5) 
39 (23.5) 

137 (82.5) 
29 (17.5) 

4 (2.4) 
1 62 (97.6) 

112 (67.5) 
54 (32.5) 

60 (36.1) 

1 06 (63.9) 

105 (63.3) 

61 (36.7) 

112 (67.5) 
54 (32.5) 

1 59 (95.8) 
7 (4.2) 

118 (71.1) 
48 (28.9) 

81 (48.8) 
85 (51.2) 

99 (59.6) 
67 (33.7) 

152 (91.6) 
14 (8.4) 



47 (68.1) 

22 (31.9) 

66 (95.7) 
3 (4.3) 

49 (71.0) 

20 (29.0) 

52 (75.4) 
1 7 (24.6) 

1 (14) 
68 (98.6) 

38 (55.1) 
31 (44.9) 

21 (30.4) 

48 (69.6) 

39 (56.5) 

30 (43.5) 

40 (58.0) 
29 (42.0) 

62 (89.9) 
7 (10.1) 

49 (71.0) 
20 (29.0) 

23 (33.3) 
46 (66.7) 

31 (44.9) 
38 (55.1) 

55 (79.7) 
14 (20.3) 



73 (35.3) 
1 34 (64.7) 

200 (96.6) 
7 (3.4) 

127 (61.3) 
80 (38.6) 

174 (84.1) 
33 (15.9) 

5 (2.4) 
202 (97.6) 

164 (79.2) 
43 (20.8) 

53 (25.6) 
1 54 (74.4) 

128 (61.8) 
79 (38.2) 

146 (70.5) 
61 (29.5) 

198 (95.7) 

9 (4.3) 

120 (58.0) 
87 (42.0) 

92 (44.4) 
115 (55.6) 

107 (51.7) 
100 (48.3) 

197 (95.2) 

10 (4.8) 



<0.001 



0.103 



0.008 



0.032 



0.027 



<0.001 



0.1 1 



0.001 



0.262 



0.262 



<0.001 



<0.001 



0.137 



<0.001 
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Table 4 Expression of metabolism-related proteins according to breast cancer phenotype* (Continued) 



LC3A in tumor 












<0.001 


Negative 


669 (90.4) 


294 (98.7) 


158 (95.2) 


68 (98.6) 


149 (72.0) 




Positive 


71 (9.6) 


4 (1.3) 


8 (4.8) 


1 (14) 


58 (28.0) 




LC3A in stroma 












<0.001 


Negative 


687 (92.8) 


267 (89.6) 


151 (91.0) 


62 (89.9) 


207 (100.0) 




Positive 


53 (7.2) 


31 (10.4) 


1 5 (9.0) 


7 (10.1) 


0 (0.0) 




LC3B in tumor 












0.013 


Negative 


475 (64.2) 


1 86 (62.4) 


1 24 (74.7) 


42 (60.9) 


1 23 (59.4) 




Positive 


265 (35.8) 


112 (37.6) 


42 (25.3) 


27 (39.1) 


84 (40.6) 




LC3B in stroma 












0.645 


Negative 


688 (93.0) 


277 (93.0) 


151 (91.0) 


65 (94.2) 


1 95 (94.2) 




Positive 


52 (7.0) 


21 (7.0) 


1 5 (9.0) 


4 (5.8) 


12 (5.8) 




Cytoplasmic p62 in tumor 












<0.001 


Negative 


274 (37.0) 


131 (44.0) 


51 (30.7) 


15 (21.7) 


77 (37.2) 




Positive 


466 (63.0) 


1 67 (56.0) 


1 1 5 (69.3) 


54 (78.3) 


130 (62.8) 




Nuclear p62 in tumor 












<0.001 


Negative 


532 (71.9) 


1 80 (60.4) 


131 (78.9) 


44 (63.8) 


177 (85.5) 




Positive 


208 (28.1) 


118(39.6) 


35 (21.1) 


25 (36.2) 


30 (14.5) 




Nuclear p62 in stroma 












0.876 


Negative 


512 (69.2) 


206 (69.1) 


1 1 5 (69.3) 


45 (65.2) 


146 (70.5) 




Positive 


228 (30.8) 


92 (30.9) 


51 (30.7) 


24 (34.8) 


61 (29.5) 





a BNIP3, BCL2/adenovirus E1B 19-kDa interacting protein 3; CAIX, carbonic anhydrase IX; Glut-1, glucose transporter 1; LC3A, microtubule-associated protein 1 light 
chain 3a; LC3B, microtubule-associated protein 1 light chain 30; MCT4, monocarboxylate transporter 4; TNBC, triple-negative breast cancer. 



Correlation between tumor metabolic phenotype and 
clinicopathologic factors 

The correlation between the metabolic phenotype of breast 
cancer and clinicopathologic parameters is summarized 
in Table 6 and Figure 1. The metabolic phenotype was 
the Warburg type (« = 298, 40.3%), the null type (n = 
326, 44.0%), the mixed type (« = 62, 8.4%) and the reverse 
Warburg type (n = 54, 7.3%). Histologic grade was highest 
in the mixed type and lowest in the null type (P < 0.001). 
The mixed type had the highest percentage of ER and PR 
negativity, and the null type had the highest percentage of 
ER and PR positivity (P < 0.001). The Warburg type had 
the highest percentage of negative HER2 status (P = 
0.006). The Warburg and mixed types comprised the 
highest percentage of TNBC, and the reverse Warburg 
and null types comprised the highest percentage of lu- 
minal A types (P < 0.001). Stromal expression of ATP 
synthase and glutaminase was high in the reverse 
Warburg and mixed types and low in the Warburg and 
null types (P < 0.001). For the status of tumor auto- 
phagy, the mixed type had the highest percentage of 
activated tumor autophagy and the null type had the 
highest percentage of nonactivated tumor autophagy 
(P < 0.001). For the status of stromal autophagy, the reverse 



Warburg and mixed types had a higher percentage of 
activation than other types (P < 0.001). Tumor expression 
of MCT4 was highest in the Warburg type and lowest in 
the null type (P < 0.001), whereas stromal expression of 
MCT4 was highest in the mixed type and lowest in the 
null type (P < 0.001). Ki-67 LI was highest in the mixed 
type and lowest in the null type (P < 0.001). 

Impact of metabolism-related proteins on patient prognosis 

The results of univariate analysis on the correlation 
between metabolism-related proteins and the clinico- 
pathologic parameters of patients are summarized in 
Table 7. Shorter disease-free survival (DFS) was associated 
with Glut-1 positivity (P = 0.010), BNIP3 negativity, tumor 
phenotype (HER2 and TNBC; P < 0.001) and tumor 
metabolic type (reverse Warburg type; P = 0.037) (Figure 2). 
Shorter overall survival (OS) was associated with Glut-1 
positivity (P = 0.023), tumor phenotype (HER2 and 
TNBC; P < 0.001) and tumor metabolic type (mixed type; 
P = 0.045) (Figure 2). Prognostic factors evaluated by 
multivariate Cox analysis included histologic grade, T stage, 
N stage, ER status, PR status, HER2 status, tumor pheno- 
type, tumor metabolic phenotype and tumor expression of 
Glut-1. The results showed that the independent factors 



Table 5 Correlations between the expression of metabolism-related proteins and clinicopathologic parameters 3 



Parameters 



Age (years, 
mean ± SD) 

Histologic grade 



ER 

Negative 
Positive 

PR 

Negative 

Positive 
HER2 

Negative 

Positive 
Tumor stage 

T1 

T2/T3 
Nodal stage 
NO 

N1/N2/N3 

Ki-67 LI (%, 
mean ± SD) 

Tumor recurrence 

Absent 

Present 
Death 

Survival 

Death 



Glut-1 in tumor 



Glut-1 in stroma 



CAIX in tumor 



CAIX in stroma 



Negative Positive 

(n = 504) (n = 236) P-value 

(%) (%) 

50.5 ± 10.7 48.1 ±11.4 0.126 

<0.001 

<0.001 

<0.001 

0.714 

<0.001 

1.932 

<0.001 
0.210 

0.420 



392(71.8) 99(41.9) 

112(22.2) 137 (58.1) 

1 28 (25.4) 1 58 (66.9) 

376(74.6) 78 (33.1) 

190(37.7) 182 (77.1) 

314(62.3) 54 (22.9) 

386(76.6) 197 (83.5) 

118(23.4) 39 (16.5) 

270 (53.6) 88 (37.3) 

234(46.4) 148(62.7) 

286(56.7) 150 (63.6) 

218(43.3) 86(36.4) 

12.7 ±14.9 29.6 ±22.1 



467 (92.7) 204 (86.4) 

37 (7.3) 32 (13.6) 

467 (92.7) 206 (87.3) 

37 (7.3) 30 (12.7) 



Negative Positive 

(n = 724) (n = 16) 

(%) (%) 



Negative Positive 



275 (38.0) 1 1 
449(62.0) 5(31.2) 



P-value (n = 520) 

(%) 



49.7 ±110.0 49.3 ±9.0 18.14 



482 (67.3) 9 (56.3) 
242 (33.4) 7 (43.7) 



0.378 



4.305 



6.741 



18.81 



16.12 



361 (49.9) 11 

363 (50.1) 5 (31.2) 

572 (79.0) 1 1 (68.8) 

152(21.0) 5(31.2) 

350 (48.3) 8 (50.0) 

374 (51.7) 8 (50.0) 

426 (58.8) 10 (62.5) 

298 (41 .2) 6 (37.5) 

18.0 ±19.3 22.5 ±14.7 7.497 



655(90.5) 16(100.0) 

69 (9.5) 0 (0.0) 

657 (90.7) 16 (100.0) 

67 (9.3) 0 (0.0) 



(n = 220) 
(%) 



8.085 



359 (69.0) 
161 (31.0) 

182 (35.0) 
338 (65.0) 

257 (49.4) 
263 (50.6) 

397 (76.3) 
123 (23.7) 

251 (48.3) 
269 (51.7) 

304 (58.5) 
216 (41.5) 
16.1 ± 1 7.5 



471 (90.6) 
49 (9.4) 

475 (91.3) 
45 (8.7) 



200 (90.9) 
20 (9.1) 

1 98 (90.0) 
22 (10.0) 



Negative Positive 

P-value (n = 627) (n = 113) 

(%) (%) 



49.7 ±11.0 49.8 ±11.0 19.67 



0.441 

1 32 (60.0) 
88 (40.0) 

0.042 

104 (47.3) 
116 (52.7) 

10.92 

115 (52.3) 

105 (47.7) 

0.294 

1 86 (84.5) 
34 (15.5) 

19.65 

107 (48.6) 
113 (51.4) 

15.62 

1 32 (60.0) 
88 (40.0) 

22.7 ±22.1 <0.001 



18.62 



11.76 



49.3 ±11.1 51.7+10.3 



427(68.1) 64(56.6) 

200(31.9) 49(43.4) 

233 (37.2) 53 (46.9) 

394 (62.8) 60 (53.1) 

311 (49.6) 61 (54.0) 

316 (50.4) 52 (46.0) 

502 (80.1) 81 (71.7) 

125 (19.9) 32 (28.3) 

296 (47.2) 62 (54.9) 

331 (52.8) 51 (45.1) 

371 (59.2) 65 (57.5) 

256 (40.8) 48 (42.5) 

17.8 ±19.7 19.5 ±16.7 



568(90.6) 103(91.2) 

59(9.4) 10(9.8) 

571 (91.1) 102 (90.3) 

56 (8.9) 1 1 (9.7) 



P-value 

0.840 
0.483 

1.113 

8.715 

1.260 

3.192 

15.87 

8.505 
17.85 

15.22 



Table 5 Correlations between the expression of metabolism-related proteins and clinicopathologic parameters 3 (Continued) 



Parameters 



ATP synthase in tumor 



Negative Positive 



ATP synthase in stroma 
Negative Positive 



Glutaminase in tumor 
Negative Positive 



Glutaminase in stroma 
Negative Positive 



(n = 30) 

(%) 



(n = 710) 

(%) 



P-value 



(n = 570) (n=170) P-value (n = 219) (n = 521) P-value (n = 495) (n = 245) P-value 

(%) (%) (%) (%) (%) (%) 



Age (years, 
mean ± SD) 

Histologic grade 



ER 

Negative 
Positive 

PR 

Negative 

Positive 
HER2 

Negative 

Positive 
Tumor stage 

T1 

T2/T3 

Nodal stage 
NO 

N1/N2/N3 

Ki-67 LI (%, 
mean ± SD) 

Tumor recurrence 

Absent 

Present 
Death 

Survival 

Death 



46.7 ± 9.3 



22 (733) 
8 (26.7) 

7 (23.3) 

23 (76.7) 

8 (26.7) 
22 (73.3) 

27 (90.0) 
3 (1 0.0) 

1 2 (40.0) 
1 8 (60.0) 

1 5 (50.0) 
1 5 (50.0) 
9.8 ± 11.7 



26 (86.7) 
4(13.3) 

27 (90.0) 
3 (1 0.0) 



49.9 ±11.0 2.478 
11.63 

469 (66.1) 
241 (33.9) 

1.827 

279 (39.3) 
431 (60.7) 

0.189 

364 (51.3) 
346 (48.7) 

3.570 

556 (78.3) 
154 (21.7) 

7.539 

346 (48.7) 
364 (51.3) 

7.266 

421 (59.3) 
289 (40.7) 
18.4 ±19.4 0.336 



645 (90.8) 
65 (9.2) 

646 (91 .0) 
64 (9.0) 



10.79 



15.68 



49.5 : 



392 (68.8) 
178 (31.2) 

210 (36.8) 
360 (63.2) 

274 (48.1) 
296 (51.9) 

472 (82.8) 
98 (17.2) 

276 (48.4) 
294 (51.6) 

342 (60.0) 
228 (40.0) 
17.3 ± 19.7 



517 (90.7) 
53 (9.3) 

520 (91.2) 
50 (8.8) 



50.4 ± 10.0 



99 (58.2) 

71 (41.8) 

76 (44.7) 
94 (55.3) 

98 (57.6) 

72 (42.4) 

111 (65.3) 
59 (34.7) 

82 (48.2) 
88 (51.8) 

94 (55.3) 
76 (44.7) 
20.8 ± 1 7.5 



1 54 (90.6) 

1 6 (9.4) 

1 53 (90.0) 

17 (10.0) 



7.539 
0.252 



1.533 



0.609 



<0.001 



21.00 



6.027 



0.798 
21.00 



13.60 



49.2 ± 1 1 .6 



147 (67.1) 
72 (32.9) 

77 (35.2) 
142 (64.8) 

1 02 (46.6) 
1 1 7 (53.4) 

1 70 (77.6) 
49 (22.4) 

100 (45.7) 
119 (54.3) 

1 35 (61 .6) 
84 (38.4) 
18.6 ±20.0 



195 (89.0) 

24 (11.0) 

1 94 (88.6) 

25 (11.4) 



49.9 ±10.7 11.90 
16.77 

344 (66.0) 
1 77 (34.0) 

4.536 

209 (40.1) 
312 (59.9) 

4.158 

270 (51.8) 
251 (48.2) 

13.08 

413 (79.3) 
108 (20.7) 

7.875 

258 (49.5) 
263 (50.5) 

7.728 

301 (57.8) 
220 (42.2) 
17.8 ±18.9 13.50 



49.5 ±11.1 50.4 ±10.7 5.124 



476 (91.4) 
45 (8.6) 

479 (91.9) 
42 (8.1) 



7.014 



3.381 



345 (69.7) 
1 50 (30.3) 

1 74 (35.2) 
321 (64.8) 

235 (47.5) 
260 (52.5) 

397 (80.2) 
98 (1 9.8) 

238 (48.1) 
257 (51.9) 



146 (59.6) 
99 (40.4) 

1 1 2 (45.7) 
133 (54.3) 

137 (55.9) 
108 (44.1) 

186 (75.9) 
59 (24.1) 

1 20 (49.0) 
125 (51.0) 



446(90.1) 225(91.8) 

49 (9.9) 20 (8.2) 

447 (90.3) 226 (92.2) 

48 (9.7) 1 9 (7.8) 



0.126 



0.735 



3.822 



18.39 



297(60.0) 139(56.7) 
198(40.0) 106(43.3) 
15.6 ± 17.2 22.9 ± 22.0 0.021 



10.56 



8.757 



Table 5 Correlations between the expression of metabolism-related proteins and clinicopathologic parameters 3 (Continued) 



Parameters 



BNIP3 in tumor 



BNIP3 in stroma 



MCT4 in tumor 



MCT4 in stroma 



Negative Positive 



(n =504) 

(%) 



(n = 236) 
(%) 



P-value 



Negative Positive 

(n = 700) (n = 40) 

(%) (%) 



Negative 

P-value (n = 540) 

(%) 



Positive 

(n = 200) 

(%) 



Negative Positive 
P-value (n = 418) (n = 322) 



(n = 418) 

(%) 



P-value 



(%) 



Age 

(years, mean ± SD) 
Histologic grade 



ER 

Negative 
Positive 

PR 

Negative 

Positive 
HER2 

Negative 

Positive 
Tumor stage 

T1 

T2/T3 

Nodal stage 
NO 

N1/N2/N3 

Ki-67 LI (%, 
mean ± SD) 

Tumor recurrence 

Absent 

Present 
Death 

Survival 

Death 



48.9 ± 10.9 



341 (67.7) 
163 (32.3) 

192 (38.1) 
312 (61.9) 

255 (50.6) 
249 (49.4) 

404 (80.2) 
100 (19.8) 

232 (46.0) 
272 (54.0) 

290 (57.5) 
214 (42.5) 
18.9 ±20.9 



445 (88.3) 
59 (11.7) 

454 (90.1) 
50 (9.9) 



51.6 ±11.0 0.042 
5.859 

1 50 (63.6) 
86 (36.4) 

14.40 

94 (39.8) 
142 (60.2) 

17.07 

1 1 7 (49.6) 
1 1 9 (50.4) 

4.410 

1 79 (75.8) 
57 (24.2) 

1.449 

1 26 (53.4) 
110 (46.6) 

6.237 

146 (61.9) 
90 (38.1) 

16.2 ±15.0 1.680 



49.5 ±10.8 53.2 ±12.7 



11.1 49.6 ±10.6 17.91 



49.6 ± 11.0 49.5 ±10.9 14.80 



226 (95.8) 
10 (4.2) 

219 (92.8) 
17(7.2) 



0.021 



5.712 



463 (66.1) 
237 (33.9) 

269 (38.4) 
431 (61.6) 

349 (49.9) 
351 (50.1) 

554 (79.1) 
146 (20.9) 

340 (48.6) 
360 (51.4) 

419 (59.9) 
281 (40.1) 
18.3 ± 19.5 



632 (90.3) 
68 (9.7) 

636 (90.9) 
64 (9.1) 



28 (70.0) 
1 2 (30.0) 

1 7 (42.5) 
23 (57.5) 

23 (57.5) 

1 7 (42.5) 

29 (72.5) 
11 (27.5) 

18 (45.0) 

22 (55.0) 

1 7 (42.5) 

23 (57.5) 
14.6 ± 13.7 



39 (97.5) 
1 (2.5) 

37 (92.5) 
3 (7.5) 



15.37 



12.99 



8.757 



6.762 



15.66 



0.630 



5.208 
3.444 



15.22 



386 (71.5) 
1 54 (28.5) 

1 74 (32.2) 
366 (67.8) 

232 (43.0) 
308 (57.0) 

422 (78.1) 
118 (21.9) 

283 (52.4) 
257 (47.6) 

317 (58.7) 
223 (41 .3) 
15.3 ± 17.8 



491 (90.9) 
49 (9.1) 

493 (91.3) 
47 (8.7) 



105 (52.5) 
95 (47.5) 

112 (56.0) 
88 (44.0) 

140 (70.0) 
60 (30.0) 

161 (80.5) 
39 (19.5) 

75 (37.5) 
125 (62.5) 



180 (90.0) 
20 (10.0) 

1 80 (90.0) 
20 (10.0) 



<0.001 



<0.001 



<0.001 



11.42 



<0.001 



18.20 



119 (59.5) 
81 (40.5) 

25.6 ±21.0 <0.001 



14.11 



11.90 



314 (75.1) 
1 04 (24.9) 

1 23 (29.4) 
295 (70.6) 

1 79 (42.8) 

239 (57.2) 

356 (85.2) 
62 (14.8) 

1 96 (46.9) 
222 (53.1) 

240 (57.4) 
1 78 (42.6) 
13.2 ± 16.5 



376 (90.0) 
42 (10.0) 

374 (89.5) 
44 (1 0.5) 



1 77 (55.0) 
145 (45.0) 

163 (50.6) 

1 59 (49.4) 

1 93 (60.0) 
1 29 (40.0) 

227 (70.5) 
95 (29.5) 

162 (50.3) 

160 (49.7) 

1 96 (60.9) 
126 (39.1) 
24.3 ± 20.7 



295 (91.6) 
27 (8.4) 

299 (92.9) 
23 (7.1) 



<0.001 



<0.001 



<0.001 



<0.001 



7.854 



7.686 



<0.001 
11.00 



2.562 



Table 5 Correlations between the expression of metabolism-related proteins and clinicopathologic parameters 3 (Continued) 



Parameters 



Cytoplasmic Beclin-1 



Nuclear Beclin-1 



LC3A in tumor 



LC3A in stroma 



Negative Positive 



(n = 406) 

(%) 



(n = 334) 

(%) 



P-value 



Negative Positive 

(n = 666) (n = 74) 

(%) (%) 



Negative Positive 



P-value 



(n = 669) 

(%) 



(n = 71) 

(%) 



Negative 

P-value (n = 687) 

(%) 



Positive 

(n = 53) 

(%) 



P-value 



Age (years, 
mean ± SD) 

Histologic grade 



ER 

Negative 
Positive 

PR 

Negative 

Positive 
HER2 

Negative 

Positive 
Tumor stage 

T1 

T2/T3 

Nodal stage 
NO 

N1/N2/N3 

Ki-67 LI (%, 
mean ± SD) 

Tumor recurrence 

Absent 

Present 
Death 

Survival 

Death 



48.6 ± 10.5 51.1 ±11.4 



272 (67.0) 
134 (33.0) 

143 (35.2) 
263 (64.8) 

205 (50.5) 
201 (49.5) 

323 (79.6) 
83 (20.4) 

175 (43.1) 
231 (56.9) 

236 (58.1) 
1 70 (41 .9) 
1 7.8 ± 1 9.4 



219 (65.6) 
1 1 5 (34.4) 

143 (42.8) 
191 (57.2) 

1 67 (50.0) 
167 (50.0) 

260 (77.8) 
74 (22.2) 

183 (54.8) 
151 (45.2) 

200 (59.9) 
134 (40.1) 
18.3 ± 19.1 



360 (88.7) 311 (93.1) 

46 (1 1 .3) 23 (6.9) 

368(90.6) 305 (91.3) 

38 (9.4) 29 (8.7) 



0.042 
14.61 



6.090 



19.76 



12.36 



0.042 



13.71 



0.042 



.1 50.8 ±9.9 8.316 



50.2 ±11.0 45.6 ± 9.8 0.021 



16.75 



427(64.1) 64(86.5) 

239 (35.9) 10 (13.5) 

261 (39.2) 25 (33.8) 

405 (60.8) 49 (66.2) 

342 (51.4) 30 (40.5) 

324 (48.6) 44 (59.5) 

533 (80.0) 50 (67.6) 

1 33 (20.0) 24 (32.4) 



318 (47.7) 
348 (52.3) 

398 (59.8) 
268 (40.2) 
19.0 ± 19.8 



600 (90.1) 
66 (9.9) 

600 (90.1) 
66 (9.9) 



40 (54.1) 
34 (45.9) 

38 (51.4) 
36 (48.6) 
9.5 ± 1 0.0 



71 (95.9) 
3 (4.1) 

73 (98.6) 
1 (1-4) 



<0.001 



8.022 



1.827 



0.336 



3.612 



<0.001 



2.877 



0.189 



470 (70.3) 21 (29.6) 

1 99 (29.7) 50 (70.4) 

226 (33.8) 60 (84.5) 

443 (66.2) 11 (15.5) 

309 (46.2) 63 (88.7) 

360 (53.8) 8 (11.3) 

515 (76.9) 68 (95.8) 

154 (23.1) 3 (4.2) 



325 (48.6) 
344 (51.4) 

390 (58.3) 
279 (41 .7) 
1 5.7 ± 17.2 



610 (91.2) 
59 (8.8) 

610 (91.2) 
59 (8.8) 



33 (46.5) 
38 (53.5) 

46 (64.8) 
25 (35.2) 
39.6 ± 23.1 



61 (85.9) 
10 (14.1) 



<0.001 



<0.001 



<0.001 



<0.001 



6.552 



<0.001 



4.053 



10.75 



49.7 ±11.0 49.6 ±9.2 



449 (65.4) 42 (79.2) 

238 (34.6) 1 1 (20.8) 

278 (40.5) 8 (15.1) 

409 (59.5) 45 (84.9) 

360 (52.4) 1 2 (22.6) 

327 (47.6) 41 (77.4) 

544 (79.2) 39 (73.6) 

143 (20.8) 14 (26.4) 



8 (11.3) 



323 (47.0) 
364 (53.0) 

406 (59.1) 
281 (40.9) 
1 8.7 ± 19.7 



622 (90.5) 
65 (9.5) 

623 (90.7) 
64 (9.3) 



35 (66.0) 
18 (34.0) 



19.65 
1.029 



<0.001 



<0.001 



8.022 



0.210 



16.23 



30 (56.6) 
23 (43.4) 

10.4 ± 9.6 0.032 



49 (92.5) 
4 (7.5) 



9.765 



50 (94.3) 
3 (5.7) 



Table 5 Correlations between the expression of metabolism-related proteins and clinicopathologic parameters 3 (Continued) 



Parameters LC3B in tumor LC3B in stroma Cytoplasmic p62 in tumor Nuclear p62 in tumor Nuclear p62 in stroma 





Negative 


Positive 




Negative 


Positive 




Negative 


Positive 




Negative 


Positive 




Negative 


Positive 






(n = 475) 

(%) 


(n = 265) 
(%) 


P-value 


(n = 688) 

(%) 


(n = 52) 
(%) 


P-value 


(n = 274) 

(%) 


(n = 466) 

(%) 


P-value 


(n = 532) 
(%) 


(n = 208) 

(%) 


P-value 


(n = 512) 

(%) 


(n = 228) 

(%) 


P-value 


Age 


49.4 ± 10.4 


50.4 ± 1 2.0 


4.599 


49.6 ±11.0 


51.1 ± 10.6 


15.77 


49.4 ± 10.3 


49.9 ± 11.4 


10.35 


49.4 ± 10.9 


50.6 ± 11.2 


3.717 


49.5 ±11.2 


50.3 ± 10.4 


6.741 


(years, mean ± SD) 
































Histologic grade 






3.528 






9.450 






0.210 






<0.001 






18.18 


l/ll 


324 (68.2) 


1 67 (63.0) 




459 (66.7) 


32 (61.5) 




1 98 (72.3) 


293 (62.9) 




322 (60.5) 


169 (81.3) 




341 (66.6) 


1 50 (65.8) 




III 


151 (31.8) 


98 (37.0) 




229 (33.3) 


20 (38.5) 




76 (27.7) 


173 (37.1) 




210 (39.5) 


39 (18.8) 




171 (33.4) 


78 (34.2) 




ER 






6.279 






7.980 






2.877 






<0.001 






14.34 


Negative 


1 73 (36.4) 


1 1 3 (42.6) 




269 (39.1) 


17 (32.7) 




96 (35.0) 


1 90 (40.8) 




227 (42.7) 


59 (28.4) 




195 (38.1) 


91 (39.9) 




Positive 


302 (63.6) 


152 (57.4) 




419 (60.9) 


35 (67.3) 




1 78 (65.0) 


276 (59.2) 




305 (57.3) 


149 (71.6) 




317 (61.9) 


137 (60.1) 




PR 






10.29 






5.292 






0.105 






<0.001 






17.05 


Negative 


234 (49.3) 


138 (52.1) 




350 (50.9) 


22 (42.3) 




119 (43.4) 


253 (54.3) 




293 (55.1) 


79 (38.0) 




259 (50.6) 


1 1 3 (49.6) 




Positive 


241 (50.7) 


1 27 (47.9) 




338 (49.1) 


30 (57.7) 




1 55 (56.6) 


213 (45.7) 




239 (44.9) 


1 29 (62.0) 




253 (49.4) 


1 1 5 (50.4) 




HER2 






12.07 






12.55 






<0.001 






16.04 






13.16 


Negative 


371 (78.1) 


212 (80.0) 




540 (78.5) 


43 (82.7) 




238 (86.9) 


345 (74.0) 




421 (79.1) 


162 (77.9) 




406 (79.3) 


1 77 (34.6) 




Positive 


104 (21.9) 


53 (20.0) 




148 (21.5) 


9(17.3) 




36 (13.1) 


121 (26.0) 




111 (20.9) 


46 (22.1) 




1 06 (20.7) 


51 (22.4) 




Tumor stage 






0.357 






14.00 






3.570 






0.189 






0.189 


T1 


214 (45.1) 


144 (54.3) 




331 (48.1) 


27 (51.9) 




142 (51.8) 


216 (46.4) 




241 (45.3) 


1 1 7 (56.2) 




231 (45.1) 


127 (55.7) 




T2/T3 


261 (54.9) 


121 (45.7) 




357 (51.9) 


25 (48.1) 




1 32 (48.2) 


250 (53.6) 




291 (54.7) 


91 (43.8) 




281 (54.9) 


101 (44.3) 




Nodal stage 






6.531 






13.90 






3.969 






5.901 






10.87 


NO 


273 (57.5) 


163 (61.5) 




407 (59.2) 


29 (55.8) 




1 70 (62.0) 


266 (57.1) 




320 (60.2) 


116 (55.8) 




306 (59.8) 


130 (57.0) 




N1/N2/N3 


202 (42.5) 


102 (38.5) 




281 (40.8) 


23 (44.2) 




104 (38.0) 


200 (42.9) 




212 (39.8) 


92 (44.2) 




206 (40.2) 


98 (43.0) 




Ki-67 LI (%, 
mean ± SD) 


1 8.2 ± 19.8 


1 7.8 ± 18.2 


16.10 


1 8.4 ± 1 9.4 


1 8.9 ± 17.6 


15.77 


1 6.0 ± 1 9.2 


1 9.3 ± 19.2 


0567 


21.4 ±20.8 


9.5 ± 10.6 


<0.001 


18.5 ± 19.7 


17.0 ±18.1 


7.161 


Tumor recurrence 






18.81 






9.765 






18.83 






3.339 






19.80 


Absent 


430 (90.5) 


241 (90.9) 




622 (90.4) 


49 (94.2) 




248 (90.5) 


423 (90.8) 




477 (89.7) 


1 94 (93.3) 




464 (90.6) 


207 (90.8) 




Present 


45 (9.5) 


24 (9.1) 




66 (9.6) 


3 (5.8) 




26 (9.5) 


43 (9.2) 




55 (10.3) 


14 (6.7) 




48 (9.4) 


21 (9.2) 




Death 






8.925 






12.89 






14.53 






10.03 






14.28 


Survival 


435 (91.6) 


238 (89.8) 




624 (90.7) 


49 (94.2) 




251 (91.6) 


422 (90.6) 




481 (90.4) 


192 (92.3) 




467 (91.2) 


206 (90.3) 




Death 


40 (8.4) 


27 (10.2) 




64 (9.3) 


3 (5.8) 




23 (8.4) 


44 (9.4) 




51 (9.6) 


16(7.7) 




45 (8.8) 


22 (9.7) 





a P-values are corrected for multiple testing using the Bonferroni correction. BNIP3, BCL2/adenovirus E1B 19-kDa interacting protein 3; CAIX, carbonic anhydrase IX; Glut-1, glucose transporter 1; LC3A, microtubule-associated 
protein 1 light chain 3a; LC3B, microtubule-associated protein 1 light chain 3(3; LI, labeling index; MCT4, monocarboxylate transporter 4; TNBC, triple-negative breast cancer. 
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Table 6 Clinicopathologic characteristics of patients according to metabolic phenotype 3 



Parameters 



Warburg type 
(n = 298) (%) 



Reverse Warburg 
type (n = 54) (%) 



Mixed type 
(n = 62) (%) 



Null type 
(n = 326) (%) 



P-value 



Age (years, mean ± SD) 
Histologic grade 



Tumor stage 
T1 

T2/T3 
Nodal stage 
NO 

N1/N2/N3 
Estrogen receptor status 

Negative 

Positive 
Progesterone receptor status 

Negative 

Positive 
HER2 status 

Negative 

Positive 
Molecular subtype 

Luminal A 

Luminal B 

HER2 

Triple-negative 
ATP synthase in tumor 

Negative 

Positive 
ATP synthase in stroma 

Negative 

Positive 
Glutaminase in tumor 

Negative 

Positive 
Glutaminase in stroma 

Negative 

Positive 
Tumor mitochondrial status 

Dysfunctional 

Functional 
Stroma mitochondrial status 

Dysfunctional 

Functional 



48.5 ± 1 1 .7 

1 69 (56.7) 
129 (43.3) 

123 (41.3) 
1 75 (58.7) 

177 (59.3) 
121 (40.6) 

152 (51.0) 
146 (49.0) 

181 (60.7) 
1 1 7 (39.3) 

252 (84.6) 
46 (15.4) 

91 (30.5) 
58 (19.5) 
22 (7.4) 

1 27 (42.6) 

8 (2.7) 
290 (97.3) 

247 (82.9) 
51 (17.1) 

84 (28.2) 
214 (71.8) 

206 (69.1) 

92 (30.9) 

94 (31.5) 
204 (68.5) 

1 3 (4.4) 
285 (95.6) 



52.0 ± 10.2 

41 (75.9) 
13 (24.1) 

29 (53.7) 
25 (46.3) 

27 (50.0) 
27 (50.0) 

1 5 (27.8) 
39 (72.2) 

22 (40.7) 
32 (59.3) 

36 (66.7) 
18 (33.3) 

22 (40.7) 
18 (33.3) 
7 (13.0) 
7 (13.0) 

1 (1.9) 
53 (98.1) 

29 (53.7) 
25 (46.3) 

13 (24.1) 
41 (75.9) 

20 (37.0) 
34 (63.0) 

20 (37.0) 
34 (63.0) 

3 (5.6) 
51 (94.4) 



51.3 : 

23 (37.0) 

39 (72.2) 

34 (54.8) 

28 (45.2) 

40 (64.5) 

22 (35.5) 

39 (62.9) 

23 (37.1) 

40 (64.5) 
22 (35.5) 

47 (75.8) 
1 5 (24.2) 

12 (19.4) 

13 (21.0) 
10 (16.1) 
27 (43.5) 

2 (3.2) 
60 (96.8) 

33 (53.2) 

29 (46.8) 

1 7 (27.4) 
45 (72.6) 

29 (46.8) 
33 (53.2) 

26 (41.9) 

36 (58.1) 

8 (12.9) 
54 (87.1) 



10.2 



50.1 ± 10.5 

258 (79.1) 
68 (20.9) 



1 72 (52.8) 
1 54 (47.2) 

1 92 (58.9) 
134 (41.1) 

80 (24.5) 
246 (75.5) 

1 29 (39.6) 
1 97 (60.4) 

248 (76.1) 
78 (23.9) 

173 (53.1) 
77 (23.6) 
30 (9.2) 
46 (14.1) 

1 9 (5.8) 
307 (94.2) 

261 (80.1) 
65 (19.9) 

105 (32.2) 
221 (67.8) 

240 (73.6) 
86 (26.4) 

96 (29.4) 
230 (70.6) 

1 6 (4.9) 
310 (95.1) 



0.052 
<0.001 



0.017 



0.457 



<0.001 



<0.001 



0.006 



<0.001 



0.178 



<0.001 



0.512 



<0.001 



0.217 



0.055 
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Table 6 Clinicopathologic characteristics of patients according to metabolic phenotype 3 (Continued) 



Tumor autophagy status 










<0.001 


Activated 


1 68 (564) 


28 (51.9) 


45 (72.6) 


117 (35.9) 




Nonactivated 


1 30 (43.6) 


26 (48.1) 


1 7 (27.4) 


209 (64.1) 




Stroma autophagy status 










<0.001 


Activated 


9 (3.0) 


11 (20.4) 


13 (21.0) 


21 (6.4) 




Nonactivated 


289 (97.0) 


43 (79.6) 


49 (79.0) 


305 (93.6) 




MCT4 in tumor 










<0.001 


Negative 


1 80 (60.4) 


38 (70.4) 


40 (64.5) 


282 (86.5) 




Positive 


1 1 8 (39.6) 


16 (29.6) 


22 (35.5) 


44 (13.5) 




MCT4 in stroma 










<0.001 


Negative 


157 (52.7) 


22 (40.7) 


20 (32.3) 


219 (67.2) 




Positive 


141 (47.3) 


32 (59.3) 


42 (67.7) 


107 (32.8) 




Ki-67 LI (%, mean ± SD) 


24.6 ± 22.5 


13.2 ± 11.1 


25.2 ± 18.6 


11.5 ± 14.1 


<0.001 


Tumor recurrence 


38 (12.8) 


6 (11.1) 


4 (6.5) 


21 (6.4) 


0.043 


Patient's death 


36 (12.1) 


5 (9.3) 


6 (9.7) 


20 (6.1) 


0.081 



a LI, labeling index; MCT4, monocarboxylate transporter 4. 



associated with shorter DFS were ER negativity (odds ratio 
(OR) = 2.7, 95% CI = 1.7 to 4.5; P < 0.001), N stage (NO vs. 
Nl/2/3, OR = 2.3, 95% CI = 1.4 to 3.8; P = 0.001) and T 
stage (Tl vs. T2/3, OR = 2.4, 95% CI = 1.3 to 4.4; P = 0.002) 
and those associated with shorter OS were ER negativity 
(OR = 3.3, 95% CI = 2.0 to 5.5; P < 0.001) and N stage 
(NO vs. Nl/2/3, OR = 2.3, 95% CI = 1.4 to 3.8; P = 0.001). 

Western blot analysis of metabolism-related proteins in 
tumor and stroma according to tumor phenotype 

Western blot analysis was performed to investigate ex- 
pression of metabolism-related proteins in tumor and 
stroma according to the tumor phenotype. The expression 
of Glut-1 and ATP synthase was higher in HER2 and 
TNBC types, and the expression was higher in tumor than 
in stroma (Figure 3). The expression of p62 was higher in 
tumor than stroma, regardless of the tumor phenotype. 

Discussion 

In the present study, we investigated the differential 
expression of metabolism-related markers according 
to the subtypes of breast cancer. Expression of gly- 
colysis markers such as Glut-1, CAIX and MCT-4 
was highest in TNBC, which is consistent with the results 
of previous studies showing higher expression of Glut-1 
and CAIX in basal-like breast cancer [19]. The active meta- 
bolic status of a tumor can be inferred from common histo- 
logic features of TNBC, such as high-grade nuclei, high- 
grade prominent necrosis and increased mitotic activity 
[20], and this was supported by the results of IHC in the 
present study. Tumor expression of Glut-1 was associ- 
ated with higher histologic grade (P < 0.001), ER nega- 
tivity {P < 0.001), higher T stage (P < 0.001), and higher 



Ki-67 LI (P < 0.001), whereas CAIX was associated with 
higher Ki-67 LI (P < 0.001). Tumor expression of MCT4 
was associated with higher histologic grade (P < 0.001), ER 
negativity (P < 0.001), PR negativity (P < 0.001), higher T 
stage (P < 0.001) and higher Ki-67 LI (P < 0.001). In 
addition, expression of Glut-1, CAIX and MCT-4 was asso- 
ciated with factors reflecting higher metabolic status. In 
contrast, tumor expression of Glut-1, CAIX and MCT-4 
was lowest in luminal A and B. Luminal type tumors tend 
to show a lower grade, lower mitotic index and less necrosis 
than the HER2 type or TNBC, suggesting nonactive meta- 
bolic status of the tumor, which was supported by the 
results of IHC. 

The expression of Glut-1, CAIX, BNIP3, MCT4, 
LC3A, LC3B and p62 was observed not only in tumor 
cells but also in stromal cells, which has not been thor- 
oughly described in previous studies. The reverse Warburg 
effect theory suggests that tumor stroma, along with 
the tumor per se, plays a role in cancer metabolism 
[2-5]. According to this theory, metabolism in stromal 
cells occurs through glycolysis due to dysfunctional 
mitochondria caused by increased autophagy, whereas 
metabolism of tumor cells occurs through OXPHOS in 
functional mitochondria. This contrasts with the conven- 
tional Warburg effect theory, which states that glycolysis 
is the major metabolic process in tumor cells. The major 
metabolic phenotypes in this study were the Warburg type 
(40.3%) and the null type (44.0%), according to the 
metabolic processes of tumor and stromal cells. We 
found that each metabolic phenotype investigated had 
different characteristics. The mixed type had higher 
histologic grade, ER negativity, PR negativity and higher 
Ki-67 LI, in contrast to the null type, which had lower 
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Figure 1 Histologic and immunohistochemical features according to metabolic phenotypes of breast cancer. The Warburg and mixed 
types show high histologic grade, estrogen receptor (ER) negativity and high Ki-67 labeling index (LI). In contrast, the reverse Warburg and null 
types show low histologic grade, ER positivity and low Ki-67 LI. ATP synthase and glutaminase were high in the reverse Warburg and mixed types 
and low in the Warburg and null types (P < 0.001). BNIP3, BCL2/adenovirus E1B 19-kDa interacting protein 3; CAIX, carbonic anhydrase IX; GLUT1, 
glucose transporter 1; H & E, hematoxylin and eosin; MCT4, monocarboxylate transporter 4. 
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histologic grade, ER positivity, PR positivity and Ki-67 LI 
(P < 0.001). As both tumor cells and stromal cells are 
glycolytic in the mixed type and nonglycolytic in the null 
type, we speculate that the mixed type is a group of 
tumors with high metabolic activity and that the null type 
consists of tumors with lower metabolic activity. The 
results of the present study show that the mixed type had 
the highest percentage of activated autophagy, whereas 
the null type had the lowest percentage, thus supporting 



this hypothesis. In addition, different molecular subtypes of 
breast cancer were classified into different metabolic types. 
TNBC constituted the highest percentage of Warburg type 
and mixed type, whereas the luminal A type constituted the 
highest percentage of reverse Warburg type and null type 
(P < 0.001). Moreover, the Warburg type and the mixed 
type were classified into groups with a higher Ki-67 LI, 
whereas the reverse Warburg type and the null type were 
classified into groups with a lower Ki-67 LI (P < 0.001). 
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Table 7 Univariate analysis of the expression of metabolism-related proteins in breast cancers and disease-free 
survival or overall survival by logrank test 3 



Immunohistochemical 
factors 



Number of patients/ 
recurrence/death 



Disease-free survival 

P-value 



Overall survival 



Mean survival 
(95% CI), months 



Mean survival 
(95% CI), months 



P-value 



Glut-1 in tumor 

Negative 

Positive 
Glut-1 in stroma 

Negative 

Positive 
CAIX in tumor 

Negative 

Positive 
CAIX in stroma 

Negative 

Positive 
ATP synthase in tumor 

Negative 

Positive 
ATP synthase in stroma 

Negative 

Positive 
Glutaminase in tumor 

Negative 

Positive 
Glutaminase in stroma 

Negative 

Positive 
BNIP3 in tumor 

Negative 

Positive 
BNIP3 in stroma 

Negative 

Positive 
MCT4 in tumor 

Negative 

Positive 
MCT4 in stroma 

Negative 

Positive 
Cytoplasmic beclin-1 

Negative 

Positive 
Nuclear beclin-1 

Negative 

Positive 



504/37/37 
236/32/30 

724/69/67 
16/0/0 

520/49/45 
220/20/22 

627/59/56 
113/10/11 

30/4/3 
710/65/64 

570/53/50 
170/16/17 

219/24/25 
521/45/42 

495/49/48 
245/20/1 9 

504/59/50 
236/10/17 

700/68/64 
40/1/3 

540/49/47 
200/20/20 

418/42/44 
322/27/23 

406/46/38 
334/23/29 

666/66/66 
74/3/1 



128 (125 to 131) 
1 1 9 (1 1 2 to 1 26) 

n/a 
n/a 



126 



125 



123 
126 

123 
127 

123 
123 

125 
116 

125 
116 

123 
128 

124 
121 

125 



1 22 to 1 30) 
1 02 to 1 1 3) 

1 22 to 1 29) 
98 to 1 08) 

90 to 114) 
1 22 to 1 29) 

121 to 129) 
1 1 2 to 1 24) 

1 1 7 to 1 28) 

1 22 to 1 30) 

1 1 9 to 1 28) 
121 to 132) 

1 1 9 to 1 27) 
1 1 9 to 1 27) 

121 to 128) 
111 to 121) 

121 to 129) 
111 to 121) 

118 to 127) 
1 24 to 1 32) 

1 1 9 to 1 28) 
1 1 8 to 1 24) 

121 to 128) 
1 06 to 1 1 5) 



0.010 



n/a 



0.740 



0.927 



0.506 



0.783 



0.323 



0.554 



0.004 



0.191 



0.550 



0.673 



0.169 



0.157 



131 (128 to 134) 
123 (118 to 128) 

n/a 
n/a 



130 
123 



129 
116 



129 
128 



129 
122 



126 
130 



128 
131 



128 
121 



129 
126 



127 
130 



129 
126 



128 
136 



127 to 132) 
1 1 7 to 1 30) 

126 to 132) 
109 to 123) 

117 to 141) 
126 to 131) 

126 to 132) 
1 1 5 to 1 28) 

120 to 131) 

127 to 133) 

125 to 131) 
1 24 to 1 33) 

124 to 131) 

126 to 135) 

126 to 131) 
1 1 2 to 1 29) 

126 to 132) 
120 to 131) 

123 to 131) 
126 to 133) 

126 to 132) 
1 23 to 1 30) 

125 to 131) 
132 to 139) 



0.023 



n/a 



0.222 



0.496 



0.936 



0.398 



0.U 



0.596 



0.426 



0.973 



0.451 



0.262 



0.566 



0.031 
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Table 7 Univariate analysis of the expression of metabolism-related proteins in breast cancers and disease-free 
survival or overall survival by logrank test 3 (Continued) 



LC3A in tumor 








0.085 






0.299 


Negative 


669/59/59 


126 (122 to ' 


1 29) 




129 (126 to 


132) 




Positive 


71/10/8 


113 (103 to • 


1 22) 




124 (115 to 


133) 




LC3A in stroma 








0.801 






0.541 


Negative 


687/65/64 


125 (122 to • 


1 29) 




128 (126 to 


131) 




Positive 


53/4/3 


65 (62 to 68) 






66 (64 to 68 


) 




LC3B in tumor 








0.990 






0.271 


Negative 


475/45/40 


125 (121 to ' 


1 30) 




130 (127 to 


133) 




Positive 


265/24/27 


118 (113 to • 


1 23) 




125 (120 to 


130) 




LC3B in stroma 








0.481 






0.565 


Negative 


688/66/64 


125 (122 to • 


1 29) 




128 (126 to 


131) 




Positive 


52/3/3 


63 (60 to 66) 






64 (62 to 66) 




Cytoplasmic p62 in tumor 








0.958 






0.528 


Negative 


274/26/23 


121 (112 to ' 


1 29) 




129 (125 to 


133) 




Positive 


466/43/44 


126 (122 to • 


1 30) 




128 (125 to 


131) 




Nuclear p62 in tumor 








0.210 






0.646 


Negative 


532/55/51 


125 (122 to ' 


1 29) 




128 (125 to 


131) 




Positive 


208/14/16 


117 (110 to ' 


1 24) 




128 (122 to 


133) 




Nuclear p62 in stroma 








0.720 






0.387 


Negative 


512/48/45 


126 (122 to ' 


1 30) 




129 (126 to 


132) 




Positive 


228/21/22 


1 04 (99 to 1 09) 




124 (118 to 


130) 




Tumor phenotype 








<0.001 






<0.001 


Luminal A 


298/15/14 


130 (126 to • 


1 33) 




134 (131 to 


137) 




Luminal B 


1 66/1 2/1 1 


129 (124 to ' 


1 34) 




130 (124 to 


135) 




HER2 


69/11/12 


111 (1 00 to ' 


121) 




119 (108 to 


130) 




TNBC 


207/31/30 


116 (109 to ' 


1 24) 




120 (114 to 


126) 




Metabolic status 








0.037 






0.045 


Warburg type 


298/38/36 


119 (112 to' 


1 26) 




124 (11 9 to 


128) 




Reverse Warburg type 


54/6/5 


90 (83 to 96) 






113 (106 to 


121) 




Mixed type 


62/4/6 


1 05 (1 00 to ' 


111) 




1 1 2 (99 to 1 


26) 




Null type 


326/21/20 


1 29 (1 26 to ' 


133) 




132 (129 to 


136) 





a BNIP3, BCL2/adenovirus E1B 19-kDa interacting protein 3; CAIX, carbonic anhydrase IX; Glut-1, glucose transporter 1; LC3A, microtubule-associated protein 1 light 
chain 3a; LC3B, microtubule-associated protein 1 light chain 30; MCT4, monocarboxylate transporter 4; n/a, not applicable; TNBC, triple-negative breast cancer. 



This result is consistent with those of a former study on 
the reverse Warburg effect in which a luminal A breast 
cancer cell line, MCF-7, was used for in vitro study [6]. 
Thus, further in vitro studies should be carried out with 
various cell lines showing different molecular subtypes. 

We identified the expression of a mitochondrial 
metabolism-related protein such as ATP synthase and 
glutaminase in the tumor and stroma in the present 
study. Notably, stromal expression of ATP synthase 
and glutaminase was high in the reverse Warburg type 
and mixed type and low in the Warburg type and null type 
(P < 0.001). We speculate that stroma showing glycolysis 



have high mitochondrial metabolic activity, as both the 
reverse Warburg and mixed types are subtypes of the 
glycolysis type of stroma by definition. Because it has 
also been reported that certain types of tumor generate 
ATP through glycolysis as well as through mitochondrial 
OXPHOS, the hypothesis that dual types of stromal me- 
tabolism via the glycolysis pathway and the mitochondrial 
pathway should be investigated further. 

The present study shows that the Warburg type and 
mixed type consisted of metabolically active and biologic- 
ally aggressive tumors, whereas the reverse Warburg 
type and null type consisted of metabolically inactive 
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Figure 2 Disease-free survival (a) and overall survival curves (b) according to the metabolic phenotypes of breast cancer. 



and biologically nonaggressive tumors. This finding 
suggests that glycolysis of tumors significantly affects 
their metabolic and biological characteristics. The asso- 
ciation of Glut-1 with shorter DFS and OS in univariate 
analysis supports this hypothesis. 

A potential limitation of this study is the use of TMA 
cores for analysis, which may not truly represent the whole 
tumor. Although it is a reasonable contention, given the 
well-known intrinsic heterogeneity of breast cancer, this 
limitation was overcome by using two 3-mm tissue cores 
because it was previously reported that TMA with two 
0.6-mm cores were representative of standard full tissue 
sections in breast cancer [21]. 

Among the breast cancer subtypes, TNBC comprised 
28% of the total cases in this study, which is higher than 
the previously reported 12% to 24%. This difference can 
be attributed to possible differences in ethnic incidence, 
as reported previously, and to the overestimation of the 
true incidence potentially by the use of TMA containing 
part of the tumor, as we defined TNBC as all negative for 
ER, PR and HER2, which is similar to known phenomena 
of the discordance of ER, PR and HER2 expression between 
samples from core biopsy and excision [22-24]. Last, 
erroneous results of ER, PR and HER2 expression may 
affect the incidence, given that a 10% of false-negative rate 
and a 5% of false-positive rate were reported in ER expres- 
sion, whereas a 4% of false-negative and false-positive cases 
were reported in HER2 [25]. Thus, cautious interpretation 
of the expression of those markers seems crucial, as misin- 
terpretation of results may lead to the misclassification of 
the molecular subtypes. 



TNBC 
T S 



LMA 



LMB HER-2 



T S T S 



T S 



Glut-1 



CAIX 



ATP5A 



GLS1 



MCT4 



LC3A 



p62 



B-actin 



HI 



Figure 3 Expression of metabolism-related proteins in tumor 
and stroma assessed by Western blot analysis according to 
tumor phenotype. CAIX, carbonic anhydrase IX; GLS1, glutaminase 
1; Glut-1, glucose transporter 1; LMA, luminal A; LMB, luminal B; 
LC3A, microtubule-associated protein 1 light chain 3a; MCT4, 
monocarboxylate transporter 4; S, stroma; T, tumor; TNBC, 
triple-negative breast cancer. 
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Conclusion 

Breast cancer is heterogeneous in its metabolic status, 
and therefore it can be classified into various metabolic 
phenotypes. Specifically, the Warburg and mixed types 
had strong associations with TNBC, whereas the reverse 
Warburg type and the null type were associated with the 
luminal type, suggesting a correlation between metabolic 
phenotype and the biology of breast cancer. 
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